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PREFACE 

This  analysis  was  conducted  for  the  Advanced  Research  Projects 
Agency  as  part  of  RAND's  continuing  study  of  the  detection  of  nuclear 
explosions.  It  examines  in  detail  the  characteristics  of  the  hydro- 
magnetic  signal  produced  by  a  high-altitude  nuclear  explosion  and 
the  generation  mechanism  for  an  assumed  model  of  the  expanding 
debris . 

This  Memorandum  should  be  of  interest  to  those  concerned  with 
geophysical  effects  due  to  nuclear  detonations  and  with  the  detection 


of  these  explosions. 


V 


SUMMARY 


This  Memorandum  advances  a  theory  for  the  generation  of  magneto- 
hydrodynamic  (MHD)  waves  by  a  high-altitude  nuclear  explosion.  These 
waves  have  been  judged  responsible  for  oscillations  of  the  geomagnetic 
field  in  the  frequency  range  of  about  .01  to  10  Hz.  For  a  simple 
model  of  the  debris  expanding  at  a  velocity  less  than  the  local  Alfven 
speed  into  the  ionosphere-exosphere,  the  dependence  of  the  MHD  signal 
strength  on  the  yield- to-mass  ratio  of  the  device  is  obtained.  The 
system  is  demonstrated  to  radiate,  and  the  MHD  radiation,  with  its 
power  spectral  density,  is  calculated.  It  is  pointed  out  that  the 
MHD  radiation  can,  in  certain  cases,  be  very  much  less  than  the  hydro- 
magnetic  yield.  In  addition,  that  portion  of  the  hydroraagnet ic  yield 
which  does  not  appear  as  radiation  is  interpreted  physically,  and  con¬ 
servation  of  energy  is  demonstrated. 

The  theory  predicts  that  the  greater  the  yield- to-mass  ratio  of 
the  device,  the  stronger  the  MHD  signal  produced,  and,  consequently, 
the  easier  it  will  be  to  detect.  Further,  the  signal  radiated  from  a 
given  device  (with  specified  hydromagnet ic  yield  and  mass)  is  shown 
to  decrease  as  the  local  Alfven  speed  at  the  burst  height  increases, 
and  an  altitude  profile  of  the  radiation  from  a  specified  device  is 
presented  together  with  some  examples  illustrating  applications  of 
the  theory. 
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I.  INTRODUCTION 

The  existence  of  low-frequency  fluctuations  in  the  geomagnetic 
field  generated  by  both  high-altitude  nuclear  explosions  and  naturally 
occurring  phenomena  has  been  of  interest  for  some  time.  These  ground- 
level  oscillations  obviously  depend  upon  the  characteristics  of  both 
the  source  of  the  disturbance  and  the  transmission  medium.  Assuming 
a  steady-state  source,  a  number  of  authors^  have  calculated  trans¬ 
mission  coefficients  in  the  frequency  range  of  about  .01  to  10  Hz.  In 
addition,  a  time-dependent  ground- level  signal  from  an  impulse  source 
has  been  computed  using  a  one- dimensional  model  for  the  system  con¬ 
sisting  of  the  earth,  ionosphere,  and  exosphere.  ^ 5 ^  These  treatments 
show  that  the  time  behavior  of  the  measured  signal  depends  largely  on 
the  resonant  frequencies  of  the  transmission  coefficients  of  the 
system.  However,  little  can  be  said  concerning  the  amplitude  of  the 
ground  signal  from  a  high-altitude  burst  until  more  is  known  about  the 
strength  of  the  generated  hydromagnetic  signal  as  a  function  of  the 
source  height,  yield,  and  mass.  Thus,  to  supplement  the  previous  work, 
which  has  been  primarily  concerned  with  propagation  mechanisms,  the 
intent  of  this  Memorandum  is  to  examine  separately  certain  features 
of  the  signal  which  are  attributable  to  the  characteristics  of  the 
source.  In  particular,  the  question  of  what  fraction  of  the  total 
hydromagnetic  yield  is  radiated  from  a  bomb  as  hydromagnetic  energy 
will  be  considered.  However,  no  attempt  will  be  made  here  to  study 
the  propagation  problem  in  any  detail. 
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For  the  purpose  of  this  analysis,  a  simple  mathematical  model  has 
been  chosen  which  should  adequately  represent  the  gross  features  of 
the  expanding  debris  after  a  high-altitude  detonation.  A  brief  physi¬ 
cal  discription  of  the  problem  is  given  in  Section  II.  In  Section  III, 
the  appropriate  wave  equations  are  derived.  For  the  idealized  propa¬ 
gation  medium  considered,  it  is  pointed  out  that  to  lowest  order  in 
the  ratio  of  wave  frequency  to  ion-cyclotron  frequency,  cu/iA  ,  only  the 
fast,  or  isotropic,  mode  is  excited.  Using  this  same  ratio  as  an 
expansion  parameter  in  a  perturbation  expansion,  the  fast  mode  can 
then  be  used  as  a  source  term  for  exciting  the  slow  or  anisotropic 
mode.  While  Argus  measurements^’^  indicate  that  the  dominant  ground 
level  signal  contained  a  considerable  anisotropic  component,  a  more 
realistic  inhomogeneous  propagation  medium  (i. e. ,  the  ionosphere- 
exosphere)  could  couple  these  anisotropic  modes  more  strongly  to  the 
waves  considered  here,  A  detailed  comparison  between  theory  and  data 
cannot  be  made  until  the  present  treatment  is  coupled  with  a  realistic 
analysis  of  the  propagation  problem. 

In  Section  IV,  the  zero-order  electromagnetic  fields  (which  are 
derived  in  Appendix  A)  are  presented  with  their  respective  Fourier 
transforms  and  the  surface  current  density. 

Section  V  presents  a  calculation  of  the  magnetohydrodynamic  (MHD) 
energy  radiated  from  the  source,  and  of  the  corresponding  power 
spectrum.  Conservation  of  energy  is  discussed,  and  the  total  work 
done  in  expanding  the  device  is  calculated  directly  from  the  stress 
tensor  in  Appendix  B.  The  MHD  radiation  efficiency  of  the  device  is 
then  defined  as  the  ratio  of  true  MHD  radiation  to  hydromagnet ic  yield, 
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and  examples  are  given  for  some  assumed  yield-to-mass  ratios  and  for 
a  model  ionosphere-exosphere.  For  very-hi.gh-altitude  bursts,  the 
ground  level  signal  will  be  predominately  affected  by  the  radiation 
fields.  It  is  pointed  out  that,  for  a  given  hydromagnet ic  yield  Y 

H 

and  debris  expansion  velocity  V,  the  radiation  increases  monoton ical ly 

with  the  ratio  V/S,  where  S  is  the  local  Alfven  speed.  For  bursts  in 

regions  where  the  Alfven  speed  is  very  large,  the  radiated  hydromag- 

netic  energy  can  be  much  less  than  Y  .  Possible  large  variations  of 

n 

hydromagnetic  radiation  with  height  are  also  demonstrated.  For 
example,  it  is  pointed  out  that  if  a  daytime  burst  with  an  Argus 
type  yield  and  mass  had  been  exploded  at  2000  km  instead  of  at  about. 
480  km,  the  hydromagnetic  radiation  would  have  been  diminished  by 
almost  two  orders  of  magnitude.  This,  of  course,  has  strong 
implications  with  regard  to  detection. 

In  Section  VI,  the  first-order  problem  is  formulated,  and  a  dis¬ 
cussion  of  the  first-order  and  total  electromagnetic  fields  is 
presented  (these  are  derived  in  Appendix  C) . 
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II.  PHYSICAL  MODEL 

In  a  nuclear  explosion,  the  liberation  of  large  amounts  of 
energy  in  a  small  period  of  time  results  in  the  bomb  materials'  being 
heated  to  several  million  degrees,  forming  a  dense  plasma.  The 
highly  conducting  gas  initially  will  expand  freely  into  the  surround¬ 
ing  medium  at  a  rate  of  typically  hundreds  of  kilometers  per  second, 
pushing  the  earth's  magnetic  field  lines  aside.  In  this  Memorandum, 
primary  interest  will  be  in  this  initial  behavior  of  the  plasma, 
which  for  typical  bursts^ 5  occurs  for  periods  on  the  order  of  a 
fraction  of  a  second.  Later,  as  the  volume  increases,  the  surface 
stresses  due  to  the  distorted  magnetic  field  increase,  while  the  in¬ 
ternal  kinetic  pressure  decreases  due  to  the  decrease  in  ion  concen¬ 
tration,  The  restoring  force  thus  increases  with  the  radius,  and 
oscillations  could  occur.  Instabilities  and  recombination  processes 
may  also  affect  the  later  expansion  stages.  In  addition,  the  surface 
stresses  due  to  the  confining  magnetic  field  are  not  spherically  sym¬ 
metric,  and  the  plasma  will  eventually  become  elongated  in  the  direc¬ 
tion  of  the  field. 

A  large  fraction  of  the  total  yield  will  be  radiated  essentially 
instantaneously  as  X-rays,  resulting  in  ionization  being  produced  in 
the  D  and  lower  E  layers  of  the  ionosphere,  which  can  greatly  effect 
propagation  in  these  regions.  In  this  Memorandum,  primary  interest 
is  in  the  generation  rather  than  the  ionospheric  propagation  of  MHD 
waves  from  bursts  above  the  E  layer,  and  X-ray  effects  will  not  be  of 
concern  in  this  regard. 
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A  number  of  authors ^ ^ ^ ^  have  considered  the  geomagnetic 
effects  produced  by  a  static  diamagnetic  sphere  as  a  model  for  the 
plasma  cloud.  The  radius  of  expansion  has  been  approximately  estim- 
mated  by  equating  the  hydromagnet ic  yield,  Y  ,  of  the  bomb  to  the 

il 

magnetic  energy  expelled  by  the  conductor.  From  a  knowledge  of  the 
mass  of  the  bomb  products,  a  mean  expansion  velocity  and,  hence,  ex¬ 
pansion  time  may  be  estimated.  As  far  as  is  known,  no  estimates  have 

been  made  of  the  fraction  of  Y  which  is  emitted  as  MHD  radiation. 

ri 

In  this  paper,  the  plasma  is  taken  to  be  a  perfectly  conducting 
sphere  which  expands  with  uniform  velocity  V  for  a  time  r  and  then 
immediately  comes  to  rest.  The  plasma  medium  interface  is  assumed  to 
be  sharp,  and  is  the  site  of  eddy  currents  which  keep  the  geomagnetic 
field  from  penetrating  the  plasma. 

The  sphere  is  considered  to  expand  at  a  velocity  much  less  than 
the  Alfven  speed,  S,  into  a  homogeneous  plasma  immersed  in  a  uniform 
magnetic  field.  While  the  first  assumption  limits  the  yield- to-mass 
ratios  for  which  the  analysis  can  be  applied,  there  are  some  realistic 
cases  for  which  this  criterion  is  well  satisfied.  The  plasma  medium  is 
considered  to  be  infinite,  cold,  and  collisionless,  with  the  relevant 
electromagnetic  properties  represented  by  a  dielectric  tensor.  This 
representation,  and  the  mathematical  formulation  of  the  pertinent 
equations,  are  discussed  in  Section  III. 
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III.  MATHEMATICAL  FORMULATION  AND  DISCUSSION 

In  the  linear  hydromagnetic  approximation,  outside  of  a  sphere  of 
a  radius  St,  where  t  is  time  measured  after  burst,  the  medium  is  quies¬ 
cent,  and  the  magnetic  field,  B  e  ,  is  constant  in  space  and  time. 

o  — z 

Over  most  of  the  medium,  the  magnetic  field  will  be  perturbed  very 
little  from  its  equilibrium  value,  and  a  linearized  theory  is  appro¬ 
priate.  Within  the  framework  of  linear  theory,  the  equations  of 
motion  for  each  plasma  component  can  easily  be  solved  to  compute  the 
dielectric  tensor  of  the  medium. 

In  order  to  solve  the  problem  self-consistently ,  it  is  necessary 

to  express  the  currents,  or  particle  velocities,  in  terms  of  the 

fields  everywhere,  near  the  surface  of  the  sphere  as  well  as  far  from 

it.  Near  the  sphere,  however,  the  time- dependent  magnetic  field 

induced  by  the  moving  conductor  cannot  be  neglected  in  comparison  to 

B  .  In  fact,  the  induced  field  must  be  of  the  order  B  at  the  surface 
o  o 

in  order  that  the  total  radial  component  of  B_  vanish  there  (see 
Eq.  (6a)).  Such  an  approach  therefore  leads  to  nonlinear  equations 
for  the  particle  velocities  near  the  sphere. 

However,  by  using  linear  theory  near  the  sphere  as  well  as  far 
from  it,  it  is  found  that  the  fields  thus  calculated  have  the  correct 
order  of  magnitude  at  the  surface  (i.  e.  ,  |b_|  ^  B  ,  |Ej  «"«  VBq),  and 

that  the  boundary  conditions  can  be  exactly  satisfied.  The  surface 
current  density,  the  particle  trajectories,  and  hence  the  currents  in 
the  medium,  must  then  be  approximately  correct  near  the  sphere. 

However,  the  greater  the  distance  from  the  surface  of  the  sphere,  the 
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more  correct  a  linearized  description  of  the  properties  of  the  medium 
becomes.  For  distances  greater  than  approximately  two  radii,  the 
functional  dependence  of  the  fields  on  distance  must  be  quite  accurate 
since  linear  theory  then  becomes  an  excellent  approximation.  There¬ 
fore,  a  calculation  based  solely  upon  linearized  theory  would  yield 
correct  results  everywhere,  at  least  within  an  order  of  magnitude. 

The  equations  governing  the  propagation  of  waves  at  frequencies 
0  ^  cju  djj  Q  (where  CL  =  ZeB  /1L  ,  the  ion-cyclotron  frequency  in  MKS 
units,  Ze  is  the  ionic  charge,  and  M.  is  the  ionic  mass)  are  now  derived. 

The  dielectric  tensor  can  be  written  for  a  Cartesian  coordinate 

_  (12) 

system  as 


where 


(JU 


K!  =  1  + 


_£i_ 


2  2’ 

Qi  -  cu 


v  _  m _ el. 

*2  0.  .2 
1  Q. 

1 


(1) 


(2) 


Here  u)  ^  and  cju^e  are  the  ion  and  electron  plasma  frequencies,  respec¬ 
tively,  a  time  dependence  of  e  has  been  assumed,  and  terms  of  the 

order  of  electron  to  ion  mass  have  been  neglected.  The  wave  equation 
for  the  electric  field  may  be  written  as 

2 

V  x  V  x  E.  =  ~  K  «  E 
c~ 


(3) 


Transforming  to  cylindrical  coordinates  with  z  axis  through  the  center 

of  the  sphere  (see  Fig.  1),  the  dielectric  tensor  K  is  unaltered,  with 

the  rows  and  columns  now  labelled  (p,cp,z).  In  addition,  for  cju/Q.  «  1 

and  Q.  «  u>  .  «  uo  ,  K„  »  K,,K„.  On  a  relative  scale,  as  K„  —  <», 
r  pi  pe’  312  3  ’ 

the  plasma  becomes  infinitely  conducting  in  the  direction  of  the  zero- 

order  magnetic  field,  and  from  Eqs.  (1)  and  (3),  E  0.  In  this 

z 

limit,  the  cp  and  p  components  of  Eq.  (3)  yield,  respectively, 


2  2 

K  to  _  K  a)  _ 

(V  x  V  x  E)  +  — - —  E  =  - —  E 


cp 


CP 


(4a) 


2-  2 
5  E  K.oT 

_ e.  ,  _1 _ ?  _ 

2  2  Ep 

az  c  M 


¥ 


(4b) 


2  2  2 

In  the  limit  ai/Q.  ->0,  L-il+p/eB  =c/S,  and  K  *->  0,  where 
l  1  moo  2 

S  is  the  Alfven  velocity  and  is  the  mass  density  of  the  plasma. 

In  this  very  low  frequency  or  hydromagnet ic  regime,  the  medium  is 

nondispers ive ,  and  Eqs.  (4a)  and  (4b)  become  decoupled  for  the  two 

familiar  MHD  modes  of  propagation.  The  excitation  of  the  slow 

mode,  E  ,  which  propagates  along  the  field  lines,  is  dependent  upon 
P 

(13  ) 

the  source  containing  a  current  in  the  z  direction.  ' 

For  the  problem  under  consideration,  the  fields  are  essentially 
those  of  a  radiating  magnetic  dipole  with  moment  aligned  antiparallel 
to  with  surface  currents  (as  well  as  the  plasma  currents)  in  the 
cp  direction.  In  this  regime,  the  slow  mode  is  therefore  not  excited 
at  all.  Treating  cu/fh  as  a  small  parameter  in  a  perturbation  ex¬ 
pansion,  the  isotropic  (fast)  mode  is  coupled  to  the  anisotropic 


Fig.  1— Cylindrical  and  spherical  coordinate  systems  with  origin 
at  the  center  of  the  conductor  and  z  axis  parallel  to  B0 
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2 

(slow)  mode  only  to  order  (cd/Q^)  ,  while  the  reverse  is  true  to  order 

(Jd/Q,  .  E*  may  therefore  be  found  for  E  =  0,  and  then  used  as  a 
l  cp  p 

source  term  for  the  anisotropic  mode.  The  medium  remains  nondisper- 
sive  through  first  order  in  (i)/Q  ,  and  both  modes  will  propagate  at 
the  Alfven  velocity.  The  zero-  and  first-order  equations  are 


(V  ‘x  V  x  E) 


V  E  =  0 

;2  q> 


and 


S2E  2  .2 

— ^  +  %E  = 


0.  2 

i  S 


2  2 

respectively,  where  it  is  assumed  that  c  /S  »  1. 

(14) 

The  boundary  conditions  at  the  surface  are 


(5a) 


(5b) 


e  •  B  =  0 
— r  — 

e  x  [E  +  V  x  B]  -  0 


(6a) 

(6b) 


where  e  is  the  unit  radial  vector  and  V  =  Ve  is  the  expansion 
— r  —  ~r 

velocity  of  the  surface  of  the  sphere.  Equation  (6a)  implies  that 
the  magnetic  field  is  always  tangential  at  the  surface.  It  will  be 
shown  from  Eq.  (6b)  that  the  tangential  electric  field  at  the  surface 
is  -  V  x  B.  when  the  sphere  is  moving,  and  zero  after  it  stops.  The 
electric  field  therefore  "jumps"  at  the  surface  at  t  =  0  when  the 
sphere  begins  to  expand  and  at  t  =  t  after  it  stops.  The  discontin¬ 
uity  in  the  zero-order  electric  field  at  any  point  in  the  non- 
dispersive  medium  will  be,  at  most,  a  jump  discontinuity,  and  hence 
the  zero-order  vector  potential  will  be  continuous. 
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For  the  given  model,  Eq.  (5a)  is  solved  for  E  (and  its  Fourier 

CD 


transform  E  )  as  a  function  of  position 

9 

culated  in  Appendix  B.  The  results  are 


in  Appendix  A,  and  E  is  cal- 

P 

presented  in  Sections  IV  and 


V. 
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IV.  THE  ZERO-ORDER  FIELDS 


The  zero-order  vector  potential  is  derived  in  Appendix  A. 
Defining  T|  =  t  -  r/S,  and  with  3  =  V/S,  the  following  electromagnetic 
fields  are  obtained: 


B  =  B  cos 
r  o 


1  - 


3' 


(1-3)  (1+23) 


¥)3  +  3(ff 


(7a) 


B  ^  -  77  B  sin 

0  2  o 


*3 


2  + 


(1-3)  (1+2 13) 


+  3  ^  +  6 


(7b) 


_  „  B  S3  sin 

I  3  o 


V  2  (1-3)2 (1+23) 


'm)  +  2  fm 


r  / 


(7c) 


for 


and 


0  <  T)  <  (1-3) 


••■•I'- Viidh 

r  r  r 


'Hz . .Q.-fi).  t~ 

3t 


(8a) 


B  ^  =  -  B  sin  0 

e  o 


i  +  -4  +  I 


2r 


3  2 

a  ,  a  a 

+  - - I  e 


3  2  I+23  \  3  2  r 


r  r 


„  II  3  _  BS  .  .  (a 2  a 

E  =  I  8  TZn'o  s  m  0  ( —r  -  —  )  e  *- 
cp  2  o  1+2  p  \  2  r 


'Hz.  (L~P)^ 

3t 


'Hz..11zHLl 


■(8b) 


(8c) 


for 


T]  >  (l-3)r  . 
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The  spatial  dependence  of  the  wave  fields  is  similar  to  that  of 


an  oscillating  magnetic  dipole.  At  distances  much  greater  than 

1-3 


a,  only  the  l/r ,  or  radiation,  terms  of  the  wave  fields  are  of 

significance.  The  total  zero-order  induced  fields  are  plotted  versus 

time  in  Figs.  2-4.  The  magnitude  of  the  induced  fields  increases  with 

retarded  time  until  the  sphere  stops  expanding,  and  then  the  fields 

decay,  with  time  constant  (3t,  to  those  of  a  static  magnetic  dipole. 

During  expansion,  the  tangential  magnetic  field  at  the  surface 
l-3N'l  r 


^7]  =  ^-±jg£- J  )  is  constant  and  equal  to 


B  1  = .  i  B  — tt+a _ 

9s  2  o  (l-B)(l+2p) 


sin 


(9a) 


and 


9s 


(3S  B 


0s 


(9b) 


For  t  >  t,  the  surface  fields  (i.  e.  ,  Tj  =  t  -  a/S)  are  given  by 


II 


0s 


—  B  sin  0 
2  o 


1  - 


1+2  g 


(10) 


II 

with  E  =0.  The  surface  current  density  i  e  ,  correct  to  terms 

cps  S  “cp3 

of  order  V^/c^,  is  equal  to  B  /u  ,  where  c  is  the  vacuum  speed  of 

0s  o 

-7 

light,  and  p,^  =  4tt  x  10  MKS  units.  In  Fig,  5,  j  is  plotted  as  a 
function  of  time.  The  surface  current  density  remains  constant  in  time 
during  expansion,  then  "jumps"  at  the  instant  the  electric  field  at  the 
surface  vanishes,  and  finally  decays  to  just  that  value  needed  to  expel 
the  geomagnetic  field  from  the  conductor. 


Btf  -  Bo  sin# 
B0  sin  6 


Fig. 3 — The  0  component  of  induced  magnetic  field  versus  v 


SB0  sin  6 


Fig . 4 — The  induced  electric  field  versus  77 


sin 


0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0  2.2  2.4 

t  (sec ) 


Fig. 5— Surface  current  density  js  versus  time  t 
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The  Fourier  decomposition  of  the  fields  is  also  of  interest  for 
problems  involving  the  transmission  of  these  signals  through  a  dis¬ 
persive  medium,  e.  g.  ,  the  ionosphere.  Then  E  ,  the  Fourier  transform 

of  E  ,  is  given  by 

9 


CO 


CO 


E  =  r  E  (r ,  0,  t)  eiaJt  dt 
cp  J  cp 

r/S 


(11) 


where  k  =  u)/S.  Carrying  out  the  integration  yields 


E  -  H(cju)  elkr  \\  -  -^1  sin 
Lr 


9 


(12) 


where 


H(«u)  =  T 


~  B  a' 
3  o 


2  ix 


2  (1+23)  (lour)  U-ipcuT  ix 


1  ,  1  /  -ix  sin  x 

+  —  e 


and 


X  =  —  uo(  1-  3)  r. 


1  3 

For  j  (jot  j  <<  1,  H((jo)  -»  ~  B^a  ,  and 


1  3/1  ik\  .  ikr 

2  Boa  {—  -  t)  sin  6  e 

r 


(13) 


The  transform  for  j co  [  «  1/t  approaches  that  of  an  instantaneously 

3 

created  dipole  with  moment  -  2nBQa  /p,  ,  and  in  this  limit  agrees  with 
the  solution  obtained  by  Karplus. 


See  Ref.  13,  p.  803. 
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It  is  of 
from  the  model 
has  estimated 


V.  ENERGY  CONSIDERATIONS  AND  RADIATION 


interest  to  determine  the  hydromagnet ic  energy  radiated 
as  a  function  of  expansion  velocity.  Leipunskii^*^ 


that 


V  ~ 


(14) 


where  Y  is  the  hydrodynamic  yield  of  the  device  and  M  is  the  mass  of 

the  ionized  bomb  products.  Assuming  V  remains  constant,  the  total 

energy  required  to  expand  the  conductor  is  calculated  in  this  section 

as  a  function  of  its  final  radius,  and  is  equated  to  Y  .  One  can  then 

H 

estimate  the  radius  a  and  expansion  time  r  for  a  given  yield.  The 
energy  radiated  into  MHD  waves  is  also  calculated,  and  by  normalizing 
to  a  given  yield  Y^.,  the  efficiency  of  the  system  as  a  radiator  of 
hydromagnet ic  waves  can  be  defined.  A  discussion  of  the  interaction 
energy  between  the  sphere  and  the  field  is  presented  here,  and  con¬ 
servation  of  energy  is  demonstrated. 

It  is  first  observed  that  it  is  sufficient  to  retain  only  the 
zero-order  fields  in  computing  the  energy  through  first  order  in 
((ju/fh).  Then,  consider  the  radial  component  of  the  Poynting  vector 

for  the  zero-order  fields  P  =  —  E  x  B,  noting  those  terms  which  do 

—  u,  —  — 

'o 

not  vanish  when  integrated  over  an  infinite  sphere,  i. e. ,  which  fall 

2  ' 

off  no  faster  than  1/r  .  These  terms  are 


T  Q  B 
Pr  =  2 

2  p 


SB- 


.  2 

sm  Q 


o  (l-p)  (1+23) 


[2m)  +  (?; j  +  a®2] 


(15a) 
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for  0  <  T|  <  (1-(3)t,  and 


(15b) 


for  T]  >  (1-|3)t.  The  total  energy  W  passing  through  a  large  sphere  of 
radius  R  is  then  given  by 

2tf  tt  |-(  1-  j3)  t  ot  -i 

WT  =  f  dep  J  de  sin  0  R2  J  dT)  PI(r=R)  +  [  dl]  Pi:[(r=R)  (16) 

0  0  L  0  (1-P)t 


It  is  instructive  to  integrate  and  discuss  separately  the  three 
terms  in  each  of  the  brackets  in  Eq.  (15)  since  the  associated  energies 
have  different  physical  interpretations.  The  last  term  in  each  of  the 
brackets  of  Eq.  (15a)  and  Eq.  (15b)  arises  from  taking  the  vector 
product  of  the  radiation  (1/r)  electric  and  magnetic  fields  and  repre¬ 
sents  true  hydromagnetic  radiation.  Carrying  out  the  integration  over 
these  terms  yields 


for  the  hydromagnetic  wave  energy,  W  . 

K. 

The  first  two  terms  in  each  of  the  brackets  in  Eqs.  (15a)  and  (15b) 

arise  from  taking  the  vector  product  of  the  electric  field  with  the 

static  magnetic  field,  B  e  .  First  performing  the  j]  integration  over 

I  II 

the  1/r  terms  shows  that  the  contributions  from  P  and  P  exactly 

I  II 

cancel.  Hence,  the  1/r  terms  in  P  and  P  represent  energy  which 
flows  from  the  sphere  during  expansion  and  then  returns  to  it  without 
ever  being  lost  from  the  system. 
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The  integration  over  the  second  term  in  the  brackets  yields  a 
finite  contribution  to  the  energy,  which  will  temporarily  be  desig 
nated  by  W  ,  and  is  equal  to 


W 


s 


4 

3 


TT  B 

o 


2  3 
a 


d 


o 


(18) 


which  is  noted  to  be  independent  of  3.  The  total  integrated  Poynting 

vector  is  thus  +  Wn ,  but  W  is  not  associated  with  hydromagnetic 
o  K  o 

radiat ion. 

The  physical  significance  of  will  be  more  readily  understood 

after  calculating  W  ,  the  apparent  change  in  energy  stored  in  the 

magnetic  field  due  to  the  presence  of  the  sphere.  The  magnetic  energy 

2  2  , 

density  is  (B  +  )/2u,  ,  and  the  value  of  B  at  t  =  «  is  equal  to 

r  8  '  o  ~ 


B  =  B 


cos  0  (1-a^/r^)  e  -  sin  0  (l+2a^/r^)  e 


1 

■ej 


(19) 


/  \  p  2  2n 

for  r  >  a,  and  zero  for  r  <  a.  Integrating  (^1/2 [h  j  B^  -  Bq  j  over 
all  space,  and  defining  the  result  to  be  yields 


W 


B 


T,  2  3 

ttB  a 


o 

3l*o 


B 

"O 


(20) 


3 

where  M  =  -  2TTBoa  /p,  is  the  dipole  moment  of  the  sphere.  Now  the 
interaction  energy  between  the  field  and  the  conductor  is  mathemati¬ 
cally  identical  to  that  between  the  field  and  a  uniformly  magnetized 
sphere  of  magnetization  m  =  M/ 
energy  of  a  magnetized  body  in  a  uniform  magnetic  field  B^  produced 


na 
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j.  However,  the  total  interaction 
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by  fixed  currents  is  well  known  to  be^J^  W  = 
the  total  magnetic  moment  of  the  body.  Hence, 
energy  is 


-  M  •  B  ,  where  M  is 
2  —  — o 

the  real  interaction 


W 


I 


B  2  3 
ttB  a 
o 

Mv 


(21) 


which  is  precisely  W  +  W  .  The  difference  in  interpretation  between 
Eqs.  (20)  and  (21)  has  been  discussed  in  detail  by  Gilinslcy  and 
Holliday.  They  have  shown  that  in  order  to  obtain  the  correct 

energy  (Eq.  (21)),  it  is  necessary  to  compute  the  interaction  between 
the  sphere  and  a  finite  source  of  magnetic  field,  and  to  take  the 
limit  of  a  uniform  field  correctly.  The  total  interaction  energy  is 
then  W  +  W ,  where  W  can  be  interpreted  as  only  a  local  change  in 
magnetic  field  energy.  In  order  to  discuss  where  W  "goes,”  it  is 

u) 

necessary  to  assume  a  model  for  the  source  of  the  earth's  magnetic 
field.  Any  model,  however,  can  be  reduced  to  linear  combinations  of 
two  cases:  a  current  "loop"  producing  constant  magnetic  flux,  and  a 
loop  with  constant  current.  In  the  first  case,  W  as  well  as  W  will 
appear  as  stored  magnetic  energy.  In  the  constant  current  case,  the 
changing  flux  through  the  loop  will  induce  an  EMF ,  and  there  will  be 
a  net  energy  change  in  the  source  (e. g. ,  a  battery)  producing  the 
current  in  the  loop.  In  this  case,  it  can  be  shown  that  in  general 


AW,  =  -  2AW  _  _ 

battery  field 
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2  3 

Therefore,  AW.  =  -  M  *  B  =  +  2ttB  a  /p,  ,  and  the  total  change 

5  battery  —  — o  o  o 

1  2  3 

in  field  energy  is  T  —  M  *  B  =  -  ttB  a  /p,  .  In  either  of  these  two 
J  2  —  — o  o  o 

cases,  the  total  interaction  energy  is  -  77  M  •  B  =  W„  +  W  . 

L  —  — O  b  d 

A  number  of  authors  have  estimated  the  expansion  radius  by 

equating  Y  to  the  magnitude  of  W  ,  the  local  field  energy  displaced 
H  13 

by  the  conducting  sphere.  This  gives  a  correct  order  of  magnitude 

2  3 

only  because  all  of  the  energy  changes  scale  as  Bq  a  /p,^;  however, 
the  physical  reasoning  is  incorrect.  Y  should  actually  be  equated 
to  the  total  work  done  in  expanding  the  sphere  to  radius  a,  which 
can  be  calculated  by  conservation  of  energy: 


yh  ■ 


WR  +  WS 


+  wB  = 


„  2  3 
ttB  a 
o 


.(.J-.+fi.).' 


b 


o  a+2pr(i-p) 


(22) 


In  the  limit  3  ->  0 ,  in  which  there  is  negligible  radiation,  Y  is 
then  given  by  the  total  interaction  energy,  and  not  jw  j„  The  work 
may  also  be  computed  by  integrating  the  Maxwell  stresses  over  the 
surface  of  the  sphere.  This  calculation  is  carried  out  in  Appendix  B, 
and  the  result  is  shown  to  agree  with  Eq.  (22). 


Assuming  that  Y  and  B  are  known,  the  variation  of  a  with  (3 
H  o 

(e. g.  ,  varying  the  mass  M)  can  be  calculated.  In  Fig.  6, 

1/3 


ttB  /u  Y  a  versus  3  is  plotted.  The  fact  that  the  expansion 
o  o  11/ 

radius  a  decreases  with  increasing  j3  implies  that  there  is  more  radi¬ 
ation,  and  work  is  thus  expended  at  a  faster  rate  at  higher  expansion 

/  2  \1/3 

velocities.  In  Fig.  7,  S^ttB  /p,  Y r  is  plotted  as  a  function  of  3- 


The  expansion  time  t  varies  as  1/3  for  small  3,  and  decreases  more 
rapidly  for  larger  3.  The  curves  are  broken  for  3  close  to  1  because 
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in  this  region  the  use  of  linearized  theory  for  the  conductivity  is 
not  justified. 

The  ratio 


(23) 


is  the  fraction  of  energy  radiated  into  hydromagnet ic  waves  for  a 
given  yield  Y  ,  which  is  the  radiation  efficiency  of  the  system  and 
is  plotted  in  Fig.  8.  For  (3  «  1,  a  very  small  fraction  of  Y  appears 

Ll 

as  radiation,  while  for  3  close  to  1,  almost  all  of  the  yield  goes 
into  wave  energy.  Therefore,  this  model  predicts  that  for  very-high- 
altitude  bursts  where  the  ground  level  signal  is  determined  predomin¬ 
ately  by  the  radiation  field,  much  larger  hydromagnetic  signals  would 
be  produced  from  bursts  with  higher  expansion  velocities  (i.e.  , 
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where  a  =  2|3/(1-|3).  Wn  (cu)  versus  uut  is  shown  in  Fig.  9  for  |3  =  .1,  .5, 

R 

2  2 

and  .9.  W  (cju)  increases  as  u)  for  small  tu,  decreases  as  1/uu  for  large 
R 

-1 

UJ,  and  peaks  at  ojtrj  4  or  at  f  r;  .6t  ,  with  the  spectra  becoming 

flatter  in  the  vicinity  of  the  peak  as  [3  -*  1.  For  ojt  &  1  (see  Eq.  (13)), 


W  (<u) 


2  62 

B  a  cu  0 
o  .  2 

— - -  sm 

8tt|j,  S 
o 


i  %ye 

8it3  B  2S3 

O 


2 


-.2 


( 1+2 B )  (1-8) 


(1+3)' 


(ju2  sin2  0  (2  5) 


It  is  of  interest  to  compute  the  total  energy  radiated  in  the  frequency 
band  from  zero  to  some  frequency  U)^(«  4/t)  for  applications  where  band- 
limited  equipment  is  utilized.  Integrating  Eq.  (24)  over  frequencies 
and  solid  angle  yields 


f  v  2  _ 

W  =  du>  W  (u>)  dQ  =  — 
f  J  J  R  9n  B  V 

-U),-  o 


-.2 


(1+2 BV (1-6) 

.  (1+3) 3  . 


=  h  (i5p)  (1+23)2(i~3)  («JfT)3 


(2  6) 


A  few  examples  are  now  presented  to  illustrate  the  applications  of  the 
above  theory. 

EXAMPLE  1 

Leipunskii  has  estimated  the  plasma  mass  M  of  Argus  to  be 

12 

~  500  kg.  Taking  YTI  =  1  KT  =  4. 2  x  10  J,  the  expansion  velocity 
H 

V  --  ^/Yh/M  =  91  km/sec.  Then,  if  an  Argus  type  burst  were  exploded 

in  the  daytime  at  an  altitude  of  480  km,  S  250  km/sec  and  3  «=t 

-5  2 

.36.  Assuming  -  5  x  10  W/m  ,  Eq.  (22)  gives  for  the  expansion 
radius  a  =  80  km.  The  expansion  time  r  =>  a/V  =  .88  sec.  The  radiation 


WR(o>) 
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Fig. 9 — Frequency  spectrum  of  hydromagnetic  radiation 
WR  M  per  solid  angle  versus  cur 
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efficiency  W  /Y  from  Fig.  7  is  .25,  and  therefore  W  -  .25  KT.  From 
R  H.  R 

Eq.  (26)  the  energy  content  of  the  spectrum  from  0  to  is  - 
1,7  x  10  ^  gu^  KT  (o)^  «  4  sec  ^). 

EXAMPLE  2 

Consider  a  bomb  of  the  same  yield  and  mass  as  in  Example  1,  but 

exploded  at  an  altitude  of  2000  km.  Use  the  parameters  S  =  4500  km/sec 
-5  2 

and  =  5  x  10  W/m  .  Then  g  =  .02,  the  expansion  radius  a  =  87  km, 

-3 

and  r  “  .96  sec.  Tire  radiation  efficiency  is  1.2  x  10  ,  the  MED  radi- 

n  _  r  o  _  "I 

ation  W  =  1.2  x  10  ,KT,  and  WL  =  .  5  •  10  go  KT  for  uv  «  4  sec 
R  r  l  t 

EXAMPLE  3 

/io\ 

Lapp'-  '  has  used  an  estimate  of  1  KT/lb,  which  implies  an  expan- 

3 

sion  velocity  V  =  3  x  10  km/sec.  Consider  an  altitude  of  2000  km, 

and  Y  =  1  KT.  Then  g  =  .66,  a  =  64  km,  t  =  .  02  sec,  W  /Y  =  ,6, 

H  R  H 

-7  3  -1 

Wn  =  .  6  KT,  and  W  =  ,6  x  10  C0£  KT  for  uu^  «  200  sec 

R  r  r  r 

Comparing  Examples  1  and  2,  the  radiated  MHD  energy  W  for  an 

R 

Argus  type  burst  is  reduced  for  the  higher  altitude  because  of  the 
smaller  ratio  of  expansion  velocity  to  local  Alfven  speed.  In  the 
first  example,  one-fourth  of  the  hydromagnet ic  yield  appears  as  radi¬ 
ation,  while  in  the  second,  almost  all  of  the  energy  is  dissipated  by 
quasi-statically  creating  the  dipole.  Increasing  the  yield-to-mass 
ratio,  as  in  Example  3,  increases  W  and  also  the  frequency  spread 
over  which  the  radiation  occurs. 

Figure  10,  for  V  =  200  km/sec  (Y  /M  ~  .005  KT/lb),  plots  the 

ri 

radiation  efficiency  W  /Y  versus  altitude,  using  a  nighttime  sunspot 

R  ri 

(3) 


minimum  altitude  profile  for  the  Alfven  speed. 


The  two  orders  of 


3  L 


Altifude  (km) 

Fig.  10 — Radiation  efficiency  WR  /YH  versus  altitude 
using  a  nighttime  sunspot  minimum  altitude 
profile  for  the  Alfven  speed 
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magnitude  fluctuation  in  W  /Y  is  due  to  the  large  variation  in  Alfven 
speed  with  altitude.  However,  it  must  be  stressed  that  the  concept  of 
radiation  efficiency  as  defined  here  is  useful  only  as  it  relates  to 
the  strength  of  the  signal  produced  by  the  source,  and  not  to  ultimate 
detection,  which  also  depends  upon  the  propagation  medium.  The  inter¬ 
pretation  is  further  restricted  to  bursts  having  expansion  radii  less 
than  the  local  scale  height  of  the  Alfven  speed  at  the  altitude  of  de¬ 
tonation.  Finally,  for  distances  from  the  burst  such  that  the  in¬ 
equality  r  »  a  is  not  satisfied,  it  is  necessary  to  consider  all  of 
the  fields  of  Eqs.  (7)  and  (8)  as  part  of  the  signal. 
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VI.  THE  FIRST-ORDER  AND  TOTAL  FIELDS 


In  this  section,  the  generation  of  the  slow  or  anisotropic  mode 

is  analyzed  for  the  idealized  model  under  consideration.  E  ,  the 
J  cp 

solution  found  to  the  zero-order  equation,  Eq.  (5a),  is  used  as  a 

source  term  for  E  in  the  first-order  equation,  Eq.  (5b).  Equation  (5b) 
P 

is  solved  for  E  in  the  time  domain  subject  to  the  appropriate  boundary 
cond it  ions. 

The  inverse  transform  of  Eq.  (5b)  is 


2  2 
BE  .  a  E 

_ a  A _ a 

2  2  2 

az  s  at 


a3E 

J. _ _C£ 


2  3 

q.s  at 

L 


(27) 


From  Eq.  (7c),  the  righthand  side  of  Eq.  (27)  is  zero  for  0  <  n  <  (1-|3)t. 
Therefore,  to  first  order  in  (<A>/fh),  is  not  excited  at  all  until  the 

sphere  has  stopped.  The  problem  is  then  reduced  to  one  involving  time- 
independent  boundary  conditions.  Equation  (27),  or  its  transform, 

Eq.  (5b),  is  merely  the  differential  equation  governing  the  motion  of 
a  vibrating  string  with  a  distributed  source  term.  Outside  of  a  cylin¬ 
der  of  radius  P  -  r  sin  0  -  a,  the  "string"  is  infinitely  long.  For 

V~  2 

.  a  -  p  ,  and 

waves  will  be  reflected  from  the  surface. 

Equation  (5b)  is  solved  in  Appendix  C,  subject  to  the  boundary  con¬ 
ditions  above  plus  the  condition  of  outgoing  waves  at  infinity,  from  the 
Green's  function.  The  solution  is 
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ilcir 

E  =  -  k  H<^- .  p  >  a  (28a) 

p  ql  p 


a)  ,  H(co) 
Cl.  p 


ikr  ik(z+a-d) 


p  <  a 


(28b) 


and 


B 

9 


0)  k  H(a>)  „ikr 

—  _  —i— i.  cos  q  e 

Cl  S  p 


p  >  a  (29a) 


co  k  H(co) 

Cl.  S  p 

r 


cos 


0  e 


ikr 


eik(z+a-d) 


p  <  a 


(29b) 


Hence,  for  p  >  a,  the  relative  phases  of  the  waves  traveling  in 
the  plus  and  minus  z  directions  vary  with  p  in  such  a  manner  that  the 
net  effect  is  to  produce  a  spherical  wave.  The  wave  amplitudes 
decrease  with  distance  from  the  z  axis,  but  not  with  z,  and  for  suf¬ 
ficiently  large  z,  |E  |  can  be  comparable  to  Je  |,  The  total  electric 

P 

field  for  p  >  a  is  then 


E  = 


H(co)i 


ikr 


-  -  ik 
r 


sin  0  e 


cok 


-kp  Q.  sin  0  — P 

T  i 


(30) 


For  kr  »  1  the  polarization  is 


iE 

p  =  — £■  = 
E 


to/  Q. 

_ i__ 

•  ^  A 
sin  u 


(31) 


2 

For  sin  0  «  co  /  0.  and  p  >  a,  the  electric  field  is  essentially  given 
by  Eq.  (28a).  For  approximately  the  same  range  of  9  the  plane  wave 
normal  modes  are  circularly  polarized.  The  amplitude  of  these  modes 

See  Ref.  12,  p.  43. 


' — 1 1  c^J 
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can  be  related  to  the  source  fields  by  writing  Eq.  (28a)  as  the  sum 

of  a  left  and  right  circularly  polarized  wave,  each  of  magnitude 

E  I. 

P  ' 

Finally,  the  magnitude  of  the  electric  fields  inside  the  cylinder 
p  =  a  (i.  e.  ,  sin  0  <  a/r)  will  be  reduced  from  those  outside  by  the 
factor 

k[r-z-  (a-d)  ]  <  k[~  -  (a-d)] 

(valid  for  r  »  a,  k(a-d)  «  1).  The  two  waves  in  the  region  p  <  a 
thus  interfere  destructively  for  the  long  wavelengths  of  interest  here. 


VII.  CONCLUSIONS 


In  this  Memorandum,  a  theory  is  proposed  for  the  generation  of 
the  waves  responsible  for  the  prompt  magneto- telluric  effects  of  high- 
altitude  nuclear  explosions.  A  theoretical  dependence  of  the  amplitude 
of  the  generated  signal  on  the  yield-to-mass  ratio  of  the  device  is 
obtained  (for  expansion  speeds  less  than  the  Alfven  velocity)  as  well 
as  the  power  spectral  density  and  the  fraction  of  hydromagnetic  energy 
which  appears  as  MHD  radiation.  In  particular,  for  a  given  hydromag¬ 
netic  yield  Y  and  Alfven  speed  S,  the  theory  predicts  that  the  smaller 
the  mass  of  the  device,  the  stronger  the  MHD  signal  produced,  and  there¬ 
fore  the  easier  it  will  be  to  detect.  The  theory  also  predicts  a  dipole 
radiation  pattern  having  maximum  intensity  perpendicular  to  the  external 
magnetic  field,  with  most  of  the  energy  appearing  in  the  isotropic  mode. 
The  experimental  data  from  Argus,  however,  show  that  the  maximum  ground 
signal  appears  in  the  direction  of  the  field  lines  passing  through  the 
burst.  In  order  to  compare  the  signals  computed  here  with  experiment, 
it  is  necessary  also  to  consider  the  filtering  effects  of  the  earth- 
ionosphere-exosphere  system.  It  has  been  pointed  out,  for  example, 
that  inhomogeneities  in  the  propagation  medium  could  couple  the  iso¬ 
tropic  to  the  anisotropic  mode  which  propagates  along  the  field  lines. 

At  least  for  plane  waves,  the  transmission  coefficient  for  the  aniso- 

(4  5) 

tropic  mode  has  been  shown  to  dominate  that  of  the  isotropic  mode.  ’ 

If  a  similar  domination  holds  for  the  spherical  modes  generated  by  this 
model,  it  might  explain  the  apparent  discrepancy.  It  is  thus  proposed 
that,  as  far  as  is  profitable,  the  effects  of  propagation  of  these  waves 
through  the  earth- ionosphere-exosphere  system  be  incorporated  into  this 
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Appendix  A 

CALCULATION  OF  THE  ZERO- ORDER  FIELDS 


Because  the  problem  is  an  initial  value  problem  with  E^  and  its 
time  derivatives  all  initially  equal  to  zero,  Eq,  (4a)  can  be  regarded 
as  a  Laplace  transformed  equation  for  E  with  transformation  variable 

■  <P 

ito.  The  unperturbed  vector  potential  at  infinity  is  A°  =  ~  B^  r  sin  0  e  . 
The  solution  of  Eq.  (5a)  representing  outgoing  waves  and  varying  as 

(19) 

sin  0  (from  continuity  of  A  at  r  =  St)  is 


E  (r,0,ico) 

9 


K  F  (  ico) 


ioi 

Sr 


iour  /  S 
e 


sin  0 


(A-  1) 


where  F  is  an  arbitrary  function  of  its  argument  and  K  is  a  constant. 
Inverting  Eq.  (A-l)  yields 


E 

9 


(r, 0, t) 


~F(t-r/S) 

2 


F  /(t-r/S)' 


Sr 


sin  0 


where  F(t)  is  in  the  inverse  transform  of  the  function  F,  and  the  prime 
denotes  a  time  derivative.  Choosing  a  Coulomb  gauge,  and  defining 
7]  =  t-r/S,  the  corresponding  vector  potential  is 

A  =  A  -e  =  —  B  sin  8 
~  cp  -cp  2  o 

where  f(f|)  =  dF/df],  and  the  constant  K  has  been  chosen  such  that 

A  (r,0  t)  B  r  sin  0  as  r  -»  ra.  The  condition  B  =0  (Eq.  (6a))  at 

cp  2  o  r 

the  surface  of  the  sphere  implies,  with  Eq .  (A- 2) ,  that  A  vanish  there, 

T 

and  the  continuity  of  A  at  r  =  St  implies  f(0)  =  f7(0)  =  0. 


_  /ILOi 


+ 


Sr 


e 


(A-2) 
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First,  consider  the  period  from  t  =  0  to  t  =  t  when  the  sphere  is 

r  '  1 -8  \ 

expanding.  At  the  surface,  r  =  Vt ,  and  T|  =  ^  J ,  where  8  =  V/S. 

Hence,  from  Eq.  (A-2),  f  must  satisfy  the  differential  equation 

f  '(71)  +  f  f(T})  -  Sr2  =  0 
at 

r  =  Ml] 

1-8 

or 

f  "CT1)  +  “  f  (Tl)  =  (A- 3) 

v  p  71  '  d-8) 


The  required  solution  is 


1  (T|) 


_3  3  3 
S  8  Tl 


(1-8)  (1+28) 


(A-4) 


and  hence  the  vector  potential  valid  for  0  <  T|  <  ( 1— 13 )  t 5  denoted  by 

a  1  - 

A  ,  LS 

9 


(A- 5) 


The  vector  potential  after  the  sphere  has  stopped  (time  j,  radiums 
a  =  3St)  may  be  written,  as  in  Eq,  (A-2), 


■r  B  sin 
2  o 


(A- 6) 
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Again,  applying  the  condition  =  0  at  r  =  a  implies  g('f])  satisfy 


g'O])  +  f  §(T])  =  a2S 


(A- 7) 


which  has  the  solution 

g ( T|)  =  a3  +  Ce  <T^T  (A- 8) 

The  time  it  takes  the  signal  that  the  sphere  has  stopped  to  reach  a 
distance  r  is  t  =  t  +  (r-a)/S,  or  ]|  =  Equating  A"*"  and  A^  at 

j|  =  (l-p)r  shows  that  the  following  two  equations  must  be  satisfied 


a  +  Ce 


i  q-B) 
1+2  B 


(A-9a) 


Ce 


-3Ba  \  B 

Both  equations  yield  C  =  e  ,  and  therefore 


)  =  3B3S3n2 

(1-3)2(1+2B) 


l-B"’ 


(A-9b) 


g  (T)) 


i  B  p 

L  bt  j 

1+2  B 

(A- 10) 


Hence , 


II  1 

A  =  -  B  r  sin  6 
2  o 


(3\  B  Ba  sin  9/2  \ 


(a- li) 
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The  zero-order  fields  in  Eqs .  (7)  and  (8)  are  then  found  from  the 
relations 


E  = 

SA 

(A- 12a) 

9 

ST] 

B  = 
r 

“  A  cot  9 
r 

(A- 12b) 

Be 

-  -  (rA) 

r  dr 

(A- 12c) 

All  other  components  of  the  field  quantities  are  zero. 
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Appendix  B 

DIRECT  CALCULATION  OF  THE  WORK  DONE  IN  EXPANDING 
THE  SPHERE  FROM  THE  STRESS  TENSOR 

Consider  a  mathematical  sphere  of  radius  R,  where  Vt  <  R  <  Vt  =  a, 
concentric  with  the  conductor.  For  a  linear  dielectric  medium,  the 
volume  forces  may  be  expressed  in  terms  of  the  fields  by 

F  =  [  f  dV  =  Ft-  dS-t^f  (eEx  B)dV  (B-l) 


where  T  is  the  Maxwell  stress  tensor  and  is  integrated  over  the  surface 

bounding  the  volume  V.  Applied  to  the  mathematical  sphere,  Eq,  (C-l) 

states  that  the  volume  force  acting  on  the  currents  within  the  sphere 

is  equal  to  the  total  stress  transferred  across  the  surface  minus  the 

rate  at  which  electromagnetic  momentum  (e  E  x  13)  is  being  lost  from 

the  enclosed  volume.  The  force  acting  on  the  conductor  for  any  r(<  a) 

is  then  given  by  Eq.  (C-l)  evaluated  at  r  =  R  and  t  =  R/V. 

To  compute  the  first  term,  consider  a  local  Cartesian  coordinate 

system  fixed  on  sphere  R  with  axes  (x,y,z)  along  the  directions 

(e  ,  e  ,  e  )  ,  respectively.  Then  E  -=»  E  ,  B  -  B  ,  and  the  stress  trans- 
— r  — 0  -tp  z  Y 

f  2  2  1 

mi t ted  across  the  surface  element  dS  is  -  e E  /2  +  B  /2u  dS  e  . 

—  !.  oj  —x 

Substituting  the  fields  from  Eqs.  (7)  and  (8)  and  integrating  over  the 
surface  yields 


x. 


T  • 


dS 


3ttBo  (1+3)  (1+g2) 


o 


e 

— r 


for  the  stress  at  r  =  R. 


See  Ref.  14,  p.  112  ff. 


(l-3)2(l+23)2 


(B-2) 
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From  Eqs„  (7)  and  (8)  the  electromagnetic  momentum  density  is 
/ 


given  by 


n  B 

G  E  x  B  =  7—  “ 
—  —  4S  p, 


.  2 


o  (l-f3)~  (l+2p) 


w  +  2w 


2  + 


31 


( 1“  (3)  ( 1+2  P) 


3  2 

(w  +  3w  +  6w) 


(B-3) 


where  w  =  j|S/r.  The  integration  over  the  volume  between  the  conductor 
and  the  sphere  R  leads  to 


R 

(e  E  x  B)dV  =  K  £r  J  dr  ]r^w 

Vt 


,3 


L  (1-3)  (1+20) 


(w4  +  5w3  +  12w  +  12) 


+  (2w  +  4) 


where 


(B-4) 


2ttB 


K  = 


^oS  (l-p)2 ( 1+2  3 ) 


Differentiating  under  the  integral  sign:’. 


‘ JL  r 

.St 


j  (e  E  x  B)dV 


=  -  KV 


Jt=R/V 

r=R 


r=R 


w- 


JlzA 


(B-5) 


where  \  j-  denotes  the  integrand  in  Eq.  (B-4).  Carrying  out  the  algebra 
yields 


St 


(e  E  x  B)dV 


2  2  2  2 
6tt  Bq*  (3  (1+|3)  R 

- - " — -9 - y  e 

t=R/V  H0(l-P)  (1+2  (3) 

r=R 


(B-6) 
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Substituting  Eqs.  (B-2)  and  (B-6)  into  Eq.  (B-l)  yields 


F 


a+6)3 

d-p)(l+2f02 


e 

— r 


(  B  -  7  ) 


Integrating  over  R  from  0  to  a  gives  Eq.  (21)  for  the  total  work  done 
in  expanding  the  sphere. 
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Appendix  C 

CALCULATION  OF  THE  FIRST-ORDER  FIELDS 

In  cylindrical  coordinates,  Eq,  (5b)  is 


+  v2,,  = 


+  k  p,  =  F(oj)  e 


?  ?  k 

ik  (  p  +z  )  2 


Sz 


2  2 
p  +z 


(p 2+z2)k 


-  ik 


=  P(u>3p,z)  (C-l) 


~  3  2 

where  E  =  u,  and  F(ao)  =  !(iu>  /Q.S  )  H(<jo).  The  Green's  function  is  the 
P  1 

solution  of 


+  k2G  =  6(z-z  0  6  (p - p  ') 
dz 


(C-2  ) 


subject  to  the  appropriate  boundary  conditions. 

SOLUTION  FOR  p  >  a 

The  solution  of  Eq.  (C-2)  for  an  infinite  "string"  yielding  only 
outgoing  waves  at  infinity  is 


A  ■  2ke"ik(z'z,)  6<p-po 


z  <  z 


(C-3a) 


_  1  ik(z-z  ')  „  . 

G2  =2lk  6  6(f-p  > 


z  >  z 


(C-3b ) 


The  solution  of  Eq.  (C-l)  for  p  >  a  is  then 


G2  (p-p  '  ,z-z  /  )P(a),p  '  ,z  /)dz  rdp  '  +  J  J  G^(p-p  ',z-z /)P(oj,  p  '  ,z  ')dz  'dp  ' 


a  z 


ik(z-z  ') 
2ik 


-  ik(z-z  *) 


P(cu,  p,z  ')dz  '  +  J  - - -  P(uu,p,z  ')dz  '  (C-4) 
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Substituting  for  P  from  Eq.  (C-l)  and  letting  x  =  z  Vp  yields 


=  ikz  ZrP  e^Pr(l+x2)^-x]  t 

2  i+x2  LikP(L+^2)% 


-ikz  r  eikp[(1+x  }  +x^  r  1  1 

+  e  I  - r -  - ^TT  '  1  dx  (C“5) 

z/p  1+x  Likp(l+x2)2  J 


Letting  x  -*  -x  in  the  second  integral  shows  that  the  second  expression 

2  k 

is  equal  to  the  firjst  with  z  ->  -z.  Substituting  y  =  (1+x  )  2-x  the 
first  expression  can  be  written  as 


ikz  r  e ' 


1+y  L ikp (1+y  ) 


-  1  dy 


where  y  =  (l+z2/p2)2  -  z/p  =  Integrating  |  e ikpy/ ( 1+y2)  dy  by 

r 


parts  reveals  that  the  above  integral  is  equal  to 
Therefore, 


Y  eikpy 

'  2~m 

ikp (1+y  ) 


F(uj)  ikr 
^  ikp  6 


(G-6) 


SOLUTION  FOR  p  <  a 

The  retarded  solution  of  Eq.  (C-2)  for  z  >  0  with  G  =  0  at 

a  (  2  2.%  . 

z  =  d  =  (a  -p  )  is 


Gyp-p'.z-o  =  jk  J 


_  1  P  -ik(z-z  )  ik(z+z /-2d) 


6(p-p/)  d  <  z  <  z  7 
(G-7a) 
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nf  '  1  r  ik(z-z  ')  ik(z+z /-2d)  1  ,  /■.  /  .  . 

G2(p-p  ,z-z  )  =  !e  "  e  p (p- p  )  z  >  z  (G- 7b) 


The  solution  of  Eq.  (G-l)  for  p  <  a  is  then 


ikz 

e _ 

■2  ik 


J 


ihk2 '  p  dz  - 


-2ikd  T  ikz  '  , 

e  e  P  dz 


-  e 


■2  ikd 


je^' 


P  dz 


,  e^2  f  ikz  ‘  „  ,  , 

+  2iir  J  °  p  dz 


(C-8) 


Substituting  for  P  from  Eq.  (G-l)  yields 


=  F(cr) 


ikz 


Z ^  ikp [ (1+x^) 2-x] 


d/p 


1+x 


r— Vr-pl 

Likpd+xz)2  J 


dx 


+  I  -e 


ik(z-2d)  P  -ikz 


+  e 

d/p  z/p 


ikp [(1+x^) 2+x] 


1+x 


E 


(G-9) 


ikp (1+x  ) 2 

The  integrals  are  essentially  the  same  as  in  Eq.  (C-5)  and  yield 


‘] 


dx  ’ 


h  = 


ikp 


\  ikz 

e ik (r-z) 

e ik(a-d ) 

ik(z-2d) 

gik(a+d) 

J  6 

-  6 

,  2 

L1  +  (t  ) 

1  +  (^)J 

b+ml 

+  e 


■  ikz 


ik(r+z) 


1  +  ('^) 
P  / 


F(oj)  reikr  _  £ik(z+a-d) 
ikp  *•  * 


(C- 10) 


The  solution  for  z  <  0  is  identical  to  Eq.  (C-10)  with  z  -*  -z. 
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